Dynamics of d-wave Vortices: Angle-Dependent Nonlinear Hall Effect 
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We study the dynamics of vortices in d-wave superconductors using a phenomenological time- 
dependent Ginzburg-Landau equation with mixing of s and d-wave components. We present nu- 
merical simulations under an external driving current J oriented with an angle <p with respect to 
the b crystal axis, calculating the vortex motion and induced electric fields for n = oo. We find an 
intrinsic Hall effect for <p ^ which depends as ~ sin(49s), and increases non-linearly with J. 



PACS numbers: 74.60. Ge, 74.20.De, 74.25. Fy 

The discovery of high-T c superconductivity has re- 
newed the interest in the vortex state of type II super- 
conductors Q). In conventional s-wave superconductors, 
the physics of the vortex state was built upon two ba- 
sic results: (i) vortices form a triangular lattice structure 
in equilibrium as shown by Abrikosov and (ii) when 
an external current is applied there is dissipation due to 
vortex motion ||. There is now a growing experimen- 
tal and theoretical H evidence that the symmetry 
of the superconducting gap is most likely to be d x 2_ y 2 
wave. Therefore, it is now an issue of central interest 
how both the static and dynamic macroscopic properties 
of vortices are affected by the d-wave superconductivity. 

The structure of a single vortex in a d-wave supercon- 
ductor was recently studied by Volovik , Soininen et al. 
H and Schopol and Maki |5). It was found in B that a 
s-wave component is nucleated near the vortex core with 
opposite winding of the phase relative to the d-wave com- 
ponent, and a four-lobe shape of the amplitude. These 
results lead to formulate a phenomenological Ginzburg- 
Landau (GL) theory with an order parameter with two 
components corresponding to the s-wave pairing and the 
d-wave pairing respectively [p|,^0|-p^| . The main ingredi- 
ent of this theory is a mixed gradient coupling of the s 
and d components. This GL theory showed that far away 
from the core the s wave winding number changes from 
— 1 to +3 resulting in four additional nodes in the s-wave 
order parameter symmetrically placed along the ±a, ±6 
symmetry axes jl(],[ll| • Using this approach, Berlinsky et 
al. JTTJ showed that, due to the fourfold anisotropy of in- 
dividual vortices, the equilibrium lattice is non-triangular 
at large fields, as seen in YBa 2 Cu30 7 -5 §. 

The question we address in this Letter is how the four- 
lobe structure of d-wave vortices affects the transport 
properties of high-T c superconductors. We consider the 
dynamics of d-wave vortices when an external current J 
is applied. The current is oriented with an angle tp with 
respect to the b crystal axis, and we have calculated the 
induced electric field E due to vortex motion [fl3"|| . This 
is represented schematically in Fig. 1(a), as well as a pos- 
sible experimental setup for the measurement [Fig. 1(b)]. 

A vortex in motion is subject to Lorentz, viscous and 
hydrodynamic or Magnus forces which result in compo- 



nents both parallel and normal to its instantaneous ve- 
locity . This leads to transport phenomena such as the 
Hall effect in the mixed state, which is still under active 
discussion P,pT|dT8|. One can distinguish three levels of 
treatment of this problem: (i) a macroscopic approach, 
like the Bardeen-Stephen model 0; (ii) an intermediate 
approach, which considers a time-dependent Ginzburg- 
Landau (TDGL) equation; and (iii) a microscopic ap- 
proach which deals with the Bogoliubov-de Gennes equa- 
tion (BdG) [ fU3|Jri| . Up to now, all these approaches have 
been applied only to s-wave vortices. Here, we will study 
the dynamics of d-wave vortices with a time-dependent 
version of the GL equation studied in p|, pl)||ll| ] . 

The GL free energy for a d x 2_ y 2 superconductor stud- 
ied in fl |lo| , |i~l|| has an order parameter with two compo- 
nents d(r) and s(r) and is given by: 

/ = a d \d\ 2 + a s \s\ 2 + (3 d \d\ 4 + (3 sd (s* 2 d 2 + s 2 d* 2 ) (1) 

+7d|nd| 2 + 7 S |ns| 2 + 

+ lv [(n bS )*{n b d) - (n a s)*(n a d) + c.c] + h 2 /8ir 

where II = (— iV — e*A/hc) is defined in the ab plane, A 
is the vector potential and h = V x A. The critical tem- 
perature for d-wave superconductivity Tfr is assumed to 
be larger than the one corresponding to s-wave pairing, 
i.e. a d < and a d < a s , and all the other constants are 
taken to be positive. In (lO) it was obtained the relation: 
Id = 7u = 7s /2. For generality, we consider these con- 
stants to be independent with the constrain 7<j7s > I 2 - 

We assume that the TDGL equations of motion for the 
order parameters d and s are simply given as JlQj] 



d=-T^- 



S = —1 -fTT, 

OS ' 



(2) 



with r the relaxation parameter and </> the scalar poten- 
tial. We consider strongly type-II materials in the k = oo 
limit. Therefore, the magnetic field is uniform and time- 
independent, and we consider the case B || c. The system 
of equations is closed after considering current conserva- 
tion: VJ = 0, with the total current density J = J s +J n . 
The supercurrent J s is given by 
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J s = j^ld [d* (lid)} + -7, [a* (Us)] (3) 

+^{[ s * (n b d) + (u bS )* d]b 

-[s* (II a d) + (n aS )* d] a} + c.c. 

while the normal current is J„ = er„.E, with tr„ an effec- 
tive normal-state conductivity. 

These TDGL equations are solved numerically using 
a finite-difference algorithm |2Cf| in a square grid along 
the x, y axis with discretization Ax = Ay = 0.2£ and 
system size 25. 6£ x 25. 6£, with £ = ^/jd/\ad\. We take 
periodic boundary conditions for the physical quantities 
h, E,J a ,|d| and \s\. One convenient gauge choice is a 
time independent A = n v Boyx, where n v is the num- 
ber of vortices in the sample, B = &o/S and $ is 
the flux quantum, and E = — V<^>. The scalar poten- 
tial <j) can be decomposed as cj)(v) = (fio( r ) — E (i) ■ r with 
4>o (r) being periodic and satisfying V 2 0o = o^V • J s 
and E (t) = cr,^ 1 [J ex t — 1/5 J d 2 r J s (t)\, since for an 
applied external current, from global current conserva- 
tion, we have J ex t = 1/5'/ d 2 r J. To ensure the period- 
icity of J s the phases of d and s must obey the boundary 
conditions: 0(r + L x ic) = 0(r) + @ x and 9(r + L y y) = 
9(r) + gx + Q y where g — n v B^L y , dt® x = L x E ox (t) and 
dt&y = LyE oy (t). The external current is taken to be 
along the discretization y-axis, and the angle tp between 
d and x [see Fig. 1(a)] is varied by "rotating" the GL equa- 
tions. (We obtain similar results by fixing the numerical 
grid to (x, y) = (a, b) and rotating 3 ext , but the effects of 
the discretization are stronger in this case). The equa- 
tion are integrated with a time step At — 0.001h/\a d \r 
averaging over 5000 steps after a transient of 10 4 steps. 

A vortex moving with a velocity Vj, induces an average 
electric field (E) = — v/, x (B) [Q. Since from current 
conservation we have that V • E = — ct,^ 1 V • J s = 4np, 
the exchange of supercurrents into normal currents oc- 
curring in the vortex core results in a dipolar charge dis- 
tribution with momentum p _L v^, such that (E) oc p. 
Under an external current J ex t there is a Hall effect when 
the current is not parallel to (E), defining a Hall angle 
tan Oh = (E±)/(Eu). The TDGL equation has been used 
in the past to model vortex dynamics in conventional su- 
perconductors (l^Jl^]. In this case, it has been shown 
by Dorsey [^5[ that a Hall effect can be induced from 
two mechanisms: (i) by a complex relaxation parameter 
r = Ti-|-ir2 with T2 ^ (modeling a hydrodynamic con- 
tribution to the Hall effect), and (ii) by a normal state 
off-diagonal conductivity a n . xy (modeling a quasiparticle 
core contribution). We consider here the TDGL equa- 
tions (^) with L2 = and a n , X y = 0, to look for addi- 
tional contributions to the Hall effect coming from the 
lack of rotational symmetry in d-wave superconductivity. 

Let us first consider the motion of a single d-wave vor- 
tex under a current J ext — 0.6J C with J c = 2e ^° d J 7 ' i 
(the critical current of the pure d-wave solution), and 



parameter values a s /\a d \ = 1.25, 7^/7^ = l,7 a /7<j = 
2, 13 sd /I3 d = 0.5, T = l,a n = (4e 2 /h)(j d /2/3 d T), and for 
different values of tp, see Figure 2. We find that when 
3 ex t is in one of the directions of maximum symmetry 
(ip = 0°,45°,90°) the vortex moves with ± 3 ex t, and 
therefore Oh = 0. This is shown in Fig.2(a-c) for ip = 0. 
We see in Fig. 2(a) that the amplitude |d| has a square- 
like shape aligned with the ab axes as seen in the static 
case pl[ . We also see that |d| is depressed in the for- 
ward direction of motion where the vortex currents add 
to Jext- The corresponding s-wave amplitude |s| shows 
the four-lobe structure corresponding to the four satellite 
vortices fll|] , but with a strong deformation due to the 
vortex motion. Also the relative position of the satellite 
vortices is displaced backwards with respect to the static 
vortex. Furthermore, there is a finite \s\ induced at long 
distances because even in the absence of vortices a cur- 
rent Jext induces a s-wave component |2lf .The motion 
induced charge p(r) — 4 ~J V • J s is shown in Fig. 2(c). 
It corresponds to a dipole p oc E a „ in the direction of 
J ex t, and therefore 9h = 0. The d-wave state gives only 
a small additional quadrupolar contribution to E in this 
case. On the contrary, when the current J ex t is not in 
a direction of symmetry there is always a finite Hall ef- 
fect, which is maximum for <p — 22.5°, the case shown in 
Fig.2(d-f). In Fig. 2(d) and (e) we see that the amplitudes 
|d| and \s\ show now a motion-induced deformation that 
is determined both by the direction of J ex t and by the 
orientation of the ab axes (i.e. amplitude depression in 
the direction perpendicular to J ex t and four-lobe struc- 
ture in the ab orientation) . This results in vortex motion 
with a velocity vj, in a direction that is not perpendicular 
to Jext- In Fig. 2(f), we see from the plot of p(r) that the 
induced dipole p is not collinear with J ex t- Therefore, 
there is a transverse component in (E) giving Oh =/= 0. 

In general, we find that both components of the elec- 
tric field, longitudinal and transversal to the current, os- 
cillate with (p. We find that this results in a Hall angle 
with a dependence of tanOn oc sin(4(^). In Fig. 3(a) we 
show the dependence tan Oh with the angle tp of the cur- 
rent for 7^/7^ = l,7 s /7d = 2,a 8 /\a d \ = 5,(3 sd /(3 d = 0.5 
and J ex t = J = 0.6J C . In conventional superconductors, 
tan## should be current independent for small J/J c - In 
Fig. 3(b) we show the current dependence of tan Oh for 
ip = 22.5°. We see that the Hall angle depends with the 
external current as tanOn oc (J/J c ) 2 . This means that 
the Hall effect vanishes in linear response theory when 
T2 = 0, a n , X y = 0, as shown by Dorsey flfj. In a rotation- 
ally symmetric system, it should also vanish for all orders 
of J I J c . Here, we find that in the d-wave case there is 
a strong contribution to the Hall effect in the non-linear 
regime. We have also added an imaginary component 
to the dissipation parameter, T2 7^ 0. We have chosen 
T2 such to give a finite tan## for tp = similar to the 
typical experimental values ]lg| ]. We see in the insets of 
Figures 3(a) and (b) that the d-wave induced Hall effect 
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is additive to the Hall angle induced by T 2 . We have 
also studied the motion of a vortex lattice for different 
values of the field 5£ 2 /$ = 0.007,0.021,0.063, and we 
find tan 9h is nearly field independent [^3| . 

A simple understanding of this angle-dependent non- 
linear Hall effect can be obtained using the generalized 
London theory of Affleck et al |22|. The supercurrent J s 
can be written in terms of the superfluid velocity v s = 
v = V6 d - (2e/hc)A as: 



(bv b - av a )(v 2 - v 2 a ) 



le 2 r 
-(bd b - kd a ){d b v b - d a v a ) \, 



(4) 



with g 2 = \a.d\/2(3d and the small parameter e = 
(|ctd|/ 'a s ){lv I 'id) 2 • The main result here is that the v s 
is not collinear 3 S , except when oriented along the direc- 
tions of symmetry of the ab crystal. In a purely dissi- 
pative dynamics the balance of forces acting on a vortex 
line is given as ryv^ = 2ert s <f"ov s x n, with 77 a dissipation 
parameter (determined by a n and Ti |I^,|l5|) and n the 
direction of B. Here v s is the supervelocity far away from 
the vortex core induced by the driving current 3 ex t = J. 
From (|]) we obtain that the external force acting on the 
vortex line is, at first order in e, 



4> ; 



J x n 



$0 n 2ee 2 
c 2e 5 V b 



J 2 )(bJ h - kJ a ) x n. 



(5) 



Therefore, the resulting "Lorentz force" is not perpen- 
dicular to J. Since the motion induced electric field is 
E = Vi x B, we obtain a Hall angle given as 



explained from the fact that the supercurrent and the su- 
pervelocity are not collinear as shown in the generalized 
London model of [Q. Therefore, this result does not 
depend on the detail of the quasiparticle physics in the 
vortex core, which is not fully described in a TDGL treat- 
ment. Instead, this effect is clearly a consequence of the 
breaking of rotational invariance in d-wave superconduc- 
tivity. Close to T c , in mean field, the depairing current 
vanishes as J c ~ i 3/2 , with t = \T C - T\/T c , and e ~ t, 
then we have that tan## ~ t~ 2 at constant J. Therefore, 
the effect should be more noticeable close to T c , where 
the nonlinearity in currents is stronger. An experimental 
setup as shown schematically in the Fig. 1(b) would al- 
low a measurement of this angle-dependent Hall effect in 
an untwinned crystal. As we have found, one will have 
9h = 9 H + 6 H ((p). Typical experiments in YBa2Cu306+ y 
|Q, which arc at p = 0, give tan 6 H « 0.05 close to 
T c . For reasonable values of the parameters of the GL 
theory jl(| we have e « 0.2, which for J / J c = 0.5 and 
ip = 22.5° Eq. (§) gives tanfl^ « 0.02, therefore a de- 
pendence with p could be observable. In principle, the 
anisotropy between the a and b axis should also be con- 
sidered, but this can lead only to dependences ~ sin 2(p, 
at most. Therefore, we predict that a measurement of a 
Hall effect with an angular dependence ~ sin 4</> would be 
a direct evidence of the d-wave character of the High-T c 
superconductivity in the transport properties of vortices. 

We acknowledge CONICET for financial support and 
we also acknowledge discussions with L. N. Bulaevskii, 
F. de la Cruz, E. Osquiguil and G. Buscaglia. 



! nil (-),, = .U ( -j- ) sill -is. 



(6) 



with A a constant of the order of unity. Even when this 
analysis is very crude, we have obtained the same re- 
sult by a generalization of the Bardeen-Stephen model 
within a perturbative treatment of the London theory of 
Affleck et al. in the small parameter e [p3| . We have 
verified the Eq.([j]) by fitting our numerical results of 
Fig. 3. In Fig. 3(a) we have fitted the sinAip dependence 
and in Fig. 3(b) we have fitted the (J/ J c ) 2 dependence for 
low currents. In Fig. 3(c) we show the dependence with 
e = {ctd/ <*s){lv /id) 2 for different values of OLd/a s and 
Iv/ld taken independently. We see that for small e this 
is an adequate parameter combination and that the Hall 
effect is linear in e. From our numerical fits we obtain 
A k, 0.28. Therefore, a measurement of the angle and 
current dependence of tan## will give a direct estimate 
of the parameters of the GL model of Eq.(||). 

In conclusion, we have shown that in a ci-wave super- 
conductor, the four-lobe structure of vortices leads to a 
non-linear angle dependent Hall effect. The effect can be 
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FIG. 1. Schematic view of (a) the vortex with the crys- 
talline axis and the frame of reference (see text) and (b) the 
experimental setup proposed to measure the angle dependent 
Hall effect. 



FIG. 2. Contour plot of the amplitudes of d-wave ((a) 
and (d)) and s-wave ((b) and (e)) components of the order 
parameter of a moving vortex for a current J ex t = 0.6J C . The 
panels (c) (f) show the contour plot of the induced electric 
charge. The parameters are |7 S = 7^ = 7u,a s = 1.25|a,j| and 
Pd = 2/3 s d- In (a), (b) and (c) the angle between J ext and the 
crystal b axis is (p = and in (d),(e) and (f) is ip = 22.5°. 



FIG. 3. Plot of the Hall angle tan 6h (a) as function of the 
angle of the current ip for J ext = 0.6J C , (b) as a function of 
(Jext/Jc) 2 for p = 22.5°. The parameters are ^7 S = jd — 7u, 
a s — 5\ctd\ and (id = 2/3 s< j- The continuous lines are the 
fittings with Eq.(6). In the insets the squares, circles and 
triangles correspond to Y2/T-1 = 0, —0.02, —0.3. In (c) tan#H 
as a function of e = {j v /~fd) 2 \ad\/ct s for J ext = 0.6J C and 
ip = 22.5°. 
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Figure 1 J. J. Vicente Alvarez et al. 
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